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Abstract 

If P C M d is a rational polytope, then i P (t) := n Z d ) is a 

quasi-polynomial in t, called the Ehrhart quasi-polynomial of P. A 
period of ip(t) is V(P), the smallest T> E Z + such that T> ■ P has 
integral vertices. Often, T>(P) is the minimum period of ip(t), but, in 
several interesting examples, the minimum period is smaller. We prove 
that, for fixed d, there is a polynomial time algorithm which, given 
a rational polytope P C M. d and an integer n, decides whether n is a 
period of ip(t). In particular, there is a polynomial time algorithm to 
decide whether ip{t) is a polynomial. We conjecture that, for fixed d, 
there is a polynomial time algorithm to compute the minimum period 
of ip(t). The tools we use are rational generating functions. 



1 Introduction 

Given a rational polytope P C M. d (that is, a bounded subset of M. d which is 
defined by a finite collection of integer linear inequalities), define the function 

i P (t) =#Pnz d ), 

where tP is P dilated by a factor of t. Also define T> = T>(P) to be the 
smallest T> 6 Z + such that T> ■ P has integral vertices. Ehrhart proved 
[Ehr62] that i P (t) is a quasi-polynomial function with a period of V. In 
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other words, there exist polynomial functions fo(t), fi(t), . . . , fn-iit), called 
the constituents of ip(t), such that 

ip(t) = fj(n) for t = j (mod V). 

Example 1.1. P = [0, |] x [0, |] C M 2 . 
Then 

i P (t) = 



(^) 2 , fort 
(*±1) 2 ; fort odd 

□ 



even 



We know that D is a period of the quasi-polynomial ip(t). What is the 
minimum period? Certainly, it must divide T>. In most cases, in fact, it is 
exactly D. In certain interesting examples, however, the minimum period is 
smaller. 

Example 1.2. Given partitions A and /i, define the Gelfand-Tsetlin polytope 
P = Pxfi C R N , as in jDLM03j (following the classic |GC50p . where N is 
defined in terms of the lengths of A and fi. 

Then =ff(P fl Z ) is the dimension of the weight \i subspace of the irre- 
ducible representation of GL n C with highest weight A. Though 

i P {t) = #(p tA , tM n z N ) 

is a polynomial, that is, it has period one (see |KR86j ). T>(P) may be made 
arbitrarily large by suitable choice of A and /x (see [DLM03 ). 

□ 



Example 1.3. More generally, given partitions A, /i, 
|/i| = \p\, define the hive polytope P = P^ C M N as in 
of ideas from [KT99J), where is defined in terms of the length of A, /i, and 



and v such that |A| + 
BucOOj (an exposition 
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Then #(P fl Z^) is the Littlewood- Richardson coefficient c\ , defined 
to be the multiplicity of V v (the highest weight representation of GL n (C) 
corresponding to v) in V\ <g> V^. Though 

<J>(*) = #tf£*nZ") 

is a polynomial (see [DW02 ), T>(P) need not be one. □ 



Example 1.4. Given any V C Z + and any s dividing V, let P be the 
pentagon with vertices (0, 0), (0, (V, (£>, 0), and (1, 

In |M W04j . it is shown that this pentagon has T>(P) = T>, but has mini- 
mum period s. □ 



These examples raise several questions: When is the minimum period of 
i p (t) less than T>(P)1 When is ip(t) a polynomial? How can we tell what 
the minimum period of ip(t) is? These questions are wide-open, though 
MW04J gives a geometric characterization of the polygons Pel 2 such that 
ip(t) is a polynomial. Here, we attack these questions from a computational 
perspective. Can we find algorithms to answer these questions "quickly?" 

Let us be more precise. We define the input size of an algorithm to be 
the number of bits needed to encode the input into binary. In particular, the 
input size of an integer a is approximately 1 + log 2 |a| (the number of digits 
needed to write a in binary). An algorithm is called polynomial time if the 
number of steps it takes is bounded by a certain polynomial in the input size. 
Proving that an algorithm is polynomial time is generally regarded as proving 
that it is "quick," at least theoretically. See Pa p94| for general background 
on algorithms and computation complexity. 

Our algorithms will take as input a polytope P. The input size of a 
polytope defined by n linear inequalities (ci,x) < bi, where q G Z d ,6j e Z, 
is approximately 

nd +y]]og 2 \cij\ + } j log 2 \b i \. 

i,j i 

We can now state the main theorem, which we will prove in Section 4. 
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Theorem 1.5. Fix d. There is a polynomial time algorithm which, given a 
rational polytope FcR d and an integer n > 0, decides whether n is a period 
of the quasi-polynomial ip(t). 

In particular, there is a polynomial time algorithm which decides whether 
ip(t) is a polynomial (that is, whether n = 1 is a period). 



It is important that we fix d in this theorem, because problems of this 
sort become intractable if d is allowed to vary. For example, the problem of 
deciding whether P contains an integer point is NP-hard if d is not fixed. 

Naively applying Theorem 11.51 yields an algorithm to find the minimum 
period of ip(t) which, unfortunately, is not polynomial time. We would have 
to factor T>(t), which would give us a set of possible n, one of which must be 
the minimal period. We will prove the following corollary in Section 4. By 
a polynomial-time reduction, of Problem A to Problem B, we mean that, if 
there was some oracle which could solve Problem B instantaneously (more 
precisely, in the amount of time it takes to output the answer to Problem 
B), then we could use that oracle to get a polynomial time algorithm for 
Problem A. In other words, Problem A is "as easy as" Problem B. 

Corollary 1.6. Fix d. There is a polynomial-time reduction of the problem 
of finding the minimum period ofip(t), where P is a d- dimensional polytope, 
to the problem of factoring a natural number D. 



Unfortunately, the problem of factoring is probably hard. It is not known 
to be polynomial time (read, not too hard) or NP-hard (read, very hard) 
and is probably somewhere in between. Nevertheless, we make the following 
conjecture. 

Conjecture 1.7. Fix d. There is a polynomial time algorithm, which, given 
a d-dimensional polytope P , computes the minimum period of ip(t). 



The tools we will use are rational generating functions. Given a set S C 
Z rf , define the generating function 

a=(ai ,...,a ( j)GS a£S 
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Sets that are very large can sometimes be written compactly as rational 
generating functions in the form 

fiS] X) = 5 a - (l-xS(l-^)---(l-x^) ' (L8) 
where x G C d , a* G Q, ft G Z d , and fey eZ d \0. 
Example 1.9. S = {0, 1, 2, ... , n}, for some n. 
Then 

/(S; x) = 1 + x + x 2 H hi" 

1 - x n+1 
1 — x 

□ 



In Section 2, we present several tools to compute and to manipulate 
rational generating functions, most of which were proved in either |BP99j or 

[BW03j . 

Given a rational polytope P C M. d , define the generating function 

F P (t, z) = f (t) + h{t)z + ■■■ + fv-i(t)z v -\ 

where the fiit) are the constituents of ip{t). In Section 3, we will prove the 
following proposition, which will be useful in the proof of Theorem 11.51 

Proposition 1.10. Fix d. There is a polynomial time algorithm which, given 
a rational polytope P, computes Fp(t, z) as a rational generating function of 
the form M.ty) . 



Finally, in Section 4, we prove Theorem 11.51 and Corollary 11.61 

2 Rational generating function tools 

In this section, we present several tools to compute and manipulate rational 
generating functions. Except for Lemma l2~71 which is proved here, they were 
proved in either |BP99| or IBWOrij . 

First we present a tool for creating rational generating functions. 
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Theorem 2.1. (Theorem 4-4 °f L BP99l) Fix d. Then there exists a poly- 
nomial time algorithm which, for any given rational polyhedron P C W 1 , 
computes f(P H Z d ; x) in the form 

f{P n Z d ; X) = ^ C i( 1 _ x a il )(l_ x a i2 )...( 1 _ x a id )' 

where G { — 1, +1}, Pi, CLij G Z d ; awei 7^ /or a// z, j. In /act, for each i, 
an, da, ■ ■ ■ , aid is a basis of Z d . 



Example 2.2. P is the interval [0,n]. 

Then P fl Z = {0,1,2,..., n}, and we have already computed /(P HZ) = 

Once we have computed some rational generating functions, we also have 
several tools to manipulate them. 

Let /(x), with x G <C d , be a rational function in the form (jl.8|) . and let 
/ l5 / 2 , . . . , /<i G Z n be integer vectors. These vectors define the monomial map 
(p:C n ^C d given by 

z= afe, • • • , Zn) i-> (z\z /2 ,...,z' d ). 

If the image of does not lie entirely in the poles of /(x), we can define the 
function g : C n — > C by 

a(z) = /(0(z)), 

which is regular at almost every point in C n . Then g(z) is /(x) specialized 
at Xj = z \ In particular, if Zj = for all i, then g(z) is /(l, 1, . . . , 1). 

Example 2.3. 5 is a finite set. 

Then/(S;1,1,...,1) = |S|. □ 

We have the following theorem, which states that, given /(x) as a short 
rational generating function, we can find g(z) quickly. 
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Theorem 2.4. (Theorem 2.6 of \B WO 3)1 ) Let us fix k, an upper bound on the 
ki in Al.fy) . Then there exists a polynomial time algorithm, which, given /(x) 
in the form hl.ty and a monomial map : C n — > <C d such that the image of 
<p does not lie entirely in the poles o//(x), computes g(z) = f(<p(z)) in the 
form 

9 ^ = ^^(l-z b ")(l-z 6i2 )---(l-z^)' 

where s < k, /3, G Q, qi, b^ G Z n , and bij ^ for all i,j. 



Now let <7i(x) and ^(x) be Laurent power series given by 
gi (x) = ^ ft mX m and </ 2 (x) = ^™ x ™ 

Then the Hadamard product g = g\ * g 2 is defined to be the power series 

rneZ d 

Example 2.5. Si, S 2 are subsets of Z d , 

5l (x) = £ x m , and <? 2 (x) = x ™ 

Then 

(<7i*<7 2 )(x) = £ x™ 

rrae5in5 2 

□ 



More generally, we may take the Hadamard product with respect to a 
proper subset of the variables, by defining 

#i(y,z) *z 02(y,z) 

as above, except with a m and /3 m functions of y. We have the following 
theorem (which is a slightly more general version of Lemma 3.4 of [B W03j . 
but the proof is the same). 
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Theorem 2.6. Fix k, d\, and d 2 . Let y e C dl ,z E C d2 , and x = (y, z). 

Then there exists a polynomial time algorithm which, given I G Z dl+d ' 2 and 
functions 

E x Pi 
CK3-7- r - and 
(l-x^i)---(l-x^) 

E x <?i 
&7i 6^ M ^ZS 
1 — x * 1 • • • ( 1 — X D ' fc ) 

such that (l,ai), (/,&«) 7^ ; computes g = gi* z g 2 (where the Laurent power 
series are convergent on a neighborhood of (e' 1 , e /2 , . . . , e ld ) ). 



Note that I in the input of the algorithm is important. For example, if 
f( x ) = ti - > then / has two possible Laurent power series expansions 

f(x) = 1 + x + x 2 + ■ ■ ■ and f(x) = —x' 1 — x~ 2 — x~ 3 — • • • 

convergent on \x\ < 1 and \x\ > 1, respectively. In this paper, however, 
the power series we examine will actually be Laurent polynomials (which are 
convergent on all of C d ), so we will not have to worry about /. 
We present one final generating function tool. 

Lemma 2.7. Fix d and k. There is a polynomial time algorithm which, 
given rational generating functions g\ (x) and g 2 (x) in the form U.8\) which 
are known to be Laurent polynomials, decides whether g\ = g 2 . 

Remark: The lemma is also true if <7i(x) and g 2 (x) are Laurent power 
series with an infinite number of terms, but there are several complications 
which will be noted in the proof. 

Proof: Let /i(x) = #i(x) — g 2 (x). We want to decide whether h = 0. 
Suppose that 

/i(x) = c « x ". 

a&L d 

and let 

h(x) = h(x)*h(x)= $>*x a . 

a&Z d 
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We can compute h in polynomial time, using Theorem 12.61 Then h = if 
and only if h = 0. Since we know that h is a polynomial, we must simply 
check whether h(l) = Y2 a <=z d c a ls zero, which we can do in polynomial time 
using Theorem 12.41 If we did not know that h is polynomial, we would have 
to be a little more careful, and here is a sketch of what to do. We can find 
bounds M such that if c a = for all a with Halloo < M, then h is identically 
zero, using, for example, ideas from Section 5.1 of [Woo04j. Then if we take 
the Hadamard product 

/ T -M _ M+l -M _ M+l -M _ M+l\ 

V 1 — Xi 1 — X 2 1 — X d J ' 

we now have something which is known to be a Laurent polynomial, and h 
is identically zero if and only if h(l) = 0. □ 



3 Computing the generating function 

Proof of Proposition ll.lUl Computing, say, fo(t) alone would be easy, 
by interpolation. Indeed, first define 

9o (s) = f (sV). 

We may find <7o(0), <7o(l), • • • ,9o(d) in polynomial time, using Theorem 12.11 
and then interpolate, as follows. Let V be the (d + 1) x (d+ 1) Vandermonde 
matrix whose i,j entry is (i — l)- 2-1 as 1 < i,j < d + 1. Then, if go( s ) = 
a + ais + a2-5 2 + ■ ■ • + we have the following equation: 



a 




'9o(0)~ 


ad 




_9o(d)_ 



Multiplying by the inverse of V, we get the coefficients of go{s), and can then 
easily recover the coefficients of fo(t). 

We cannot, however, do this for each fi(t), sequentially, in polynomial 
time: there are T> of them, and T> may be exponential in the input size. 
Instead, we perform all T> interpolations simultaneously, using generating 
functions. 
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For < % < V-l, let 

9i {s) = f {sV + i) . 

For < j < d, let 

hj{z) = g (j) + gi (j)z + g 2 {j)z 2 + ■■■ + gv-x{j)^~ l - 
For < i < T> — 1 and < k < d, let be such that 

gi(s) = a i0 + cms + a i2 s 2 H h a id s d , 

and let 

adz) = a ok + ai k z + a 2 kZ 2 H h a V -i tk z D ~ 1 . 



Then we have that 



V ■ 



a (z) 




h (z) 


a d (z) 




h d {z) 



Therefore, if we can compute each hj(z) in polynomial time as short ratio- 
nal generating functions, then we could compute the at{z) as short rational 
generating functions by multiplying by the inverse of V. 
We compute 

hj(z) = g (j) + gi (j)z + g 2 (j)z 2 + ■■■ + gv-i(j)z v - 1 

= fo(jV) + h{jV + l)z + f 2 (jV + 2) + ■ ■ ■ + fv-i{jV + V - \)z v - x 
= ip(jV) + i P (jV + + i P (jV + 2)z 2 + ■■■ + i P (jV + V- l)z v -\ 

as follows. Given j, define the polyhedron 

Qj = {(z,y) : < z < V - 1 and y e (jV + z)P). 

Then 

f(Qy,z,y)= J2 z " E 

0<a<O-l b£(jV+a)P 

and 

h j( z ) = f{Qj\ z , !)• 

We may compute f{Qj',z,y) in polynomial time, using Theorem 12.11 and 
then perform the substitution y = 1, using Theorem 12.41 
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We have shown that we can construct the generating functions a k (z), for 
1 < k < d, in polynomial time. We must now use these generating functions 
to compute 

F P (t, z) = hit) + h{t)z + ■■■ + fv-i{t)z v - 1 . 
Since, for < j < V — 1, 

gj(s) = a j0 + ajis H h a jd s d 



and 



we have that 



fj(t) = 9j 



t-J 
V 



t — j 

fj(t) = a j0 + aj!—^- + --- + a jd 



t-J 
V 



and 



F P (t, z) = < 



CLqq + 

+ a 10 z + 



+ ... + 
_l_ ... + 



+ a-D-iflZ + ai>-i,i — ^ — z + 



+ av-i,d { — © — J z 



For < k < d, define 



b k (t,z) = a ok y— ) f a lk 
Then 



t - 1 
V 



z + ■ ■ ■ + a-v-ik 



t-V + l 
V 



JD-1 



F P (t, z) = b (t, z) + b 1 (t,z) + --- + b d (t, z). 
For each k, we will compute b k (t, z) from 

a-k\Z) = dQk + ^lfc^ + (l2k z + • • • + 0-V-l,k z 

In fact 



b k (t,z) = a k {z) -k z 



5' 



t- 1 



2+ ••• + 
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and 

X t V . ft- l\ k . ft-V + l\ k n _, 



vi + z + --- + {—d— ] z 

can be computed as a short rational generating function in polynomial time, 
by expanding all of the terms and repeatedly using the fact that, for any k, 

X/Si**' 2 ' ^ s (l^z)' Therefore we can compute the bk(t,z) and hence 

Fp(t, z) in polynomial time. □ 



4 Deciding whether n is a period 

Proof of Theorem 11.51 Given n and P, we want to decide whether n 
is a period of the quasi-polynomial ip(t). Using Proposition II . X()| we may 
compute the generating function 

F P (t, z) = f (t) + h(t)z + ■■■ + fv-i(t)z v -\ 

Define the generating function 

G n , P (t, z) = f n (t)+f n+1 (t)z+- ■ ■+fv-i(t)z v - n - 1 +f (t)z v - n +f 1 (t)z I> - n+1 +- ■ ■+f n - 1 (t)z v -\ 

Then n is a period of ip(t) if and only if Fp(t, z) = G ntP (t, z). We must show 
how to compute G„ 5 p in polynomial time. Note that 



F P (t,z) * z 



z n _ Z V 
1 - Z 



and 



F P (t,z) -k z 



1 - z 



F P {t,z) * z (z n + z n+1 + ■ • • + z v - x ) 

fn{t)z n + f n+1 (t)z n+1 + ■■■ + f V -l{t)z V - 1 

= F P (t,z) * a {l + z + --- + z n - 1 ) 
= h{t) + f l {t)z + --- + f n _ 1 {t)z n -\ 



Then 

G n , P (t,z) 



F P (t,z) -k z 



z n _ Z V 



z- n + 



F P (t, z) -k z 



JD-n 
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This can be computed in polynomial time, using Theorem 12.61 

We can decide whether F P (t, z) = G ntP (t, z) using Lemma l2~71 in polyno- 
mial time, and the proof follows. □ 

Proof of Corollary 11.61 Compute V = T>(P) by taking the least 
common multiple of the denominators of all of the coordinates of the vertices 
of P. Assume that we can find the prime factorization of T> using an oracle. 
Initialize the following loop with hq := T>. 

1. After the j'th iteration of the loop, rij is known to be a period of ip(t). 

2. For each prime factor p of rij, decide whether y is a period of ip(t). 

• If none are periods, then rij is the minimum period of ip(t), and 
we are done. 

• if — is a period of ip(t) for some p, then repeat the process with 

rij 
= J- 

This loop must terminate, because eventually we would have rij = 1. □ 
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